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ABSTRACT 



Aims. We investigate the performance of an analytic model of the 3D matter distribution, which combines perturbation 
theory with halo models, for weak-lensing statistics. 

Methods. We compare our predictions for the weak-lensing convergence power spectrum and bispectrum with numerical 
simulations and Htting formulas proposed in previous works. 

Results. We find that this model provides better agreement with simulations than published fitting formulas. This shows 
that building on systematic and physically motivated models is a promising approach. Moreover, this makes explicit the 
link between the weak-lensing statistics and the underlying properties of the 3D matter distribution, as a function of 
scale £. Thus, we obtain the contributions to the lensing power spectrum and bispectrum that arise from perturbative 
terms (complete up to one-loop) and nonperturbative terms (e.g., "1-halo" term). Finally, we show that this approach 
recovers the dependence on cosmology (for realistic scenarios). 

Key words, weak gravitational lensing; cosmology: theory - large-scale structure of Universe 



1. Introduction 

Weak gravitational lensing of background galaxies by fore- 
ground large-scale structures, the so-called cosmic shear, 
is one of the best tools to probe the nature of the 
main components of the Universe, such as dark mat- 
ter and dark energy. Thus, weak lensing has the high- 
est potential to constrain the properties of dark energy 
among other cosmological observati ons, if the systemati c 
errors are well kept under control ()Albrecht et al.l [20061) . 
To address questions about the nature of dark energy 
and the properties of gravity on cosmological scales, var- 
ious surveys are planned, su ch as the Hyper Su prime- 
Cam Weak Lensing Survey (|Mivazaki et alJ l2006[ fl. the 
Dark Energy Survey (DES fl the Large Synop tic Survey 
Telescope (LSSTfl Euclid (|Refregier et al.ll20TolPl and the 
Wide-Field Infrared Survey Telescope (WFIRStB- 

To exploit the full potential of future weak-lensing sur- 
veys, it will be important to analyze data with adequate sta- 
tistical measures and tools. Particularly, one needs to prop- 
erly take into account the correlations of the observables 
between different an gular scales and red s hifts, i.e., their 



covariance matrices ( Cooray &: Hul 120011 : iTakada fc JainI 
l2009HSato et al.ll2009Ll2011bD . Furthermore, one has to use 
an appropr iate likelihood function with given marginal dis- 
tributions (|Sato et al.ll201d■l2011a^ ■ 



^ http;/ /www. naoj.org/Projects/HSC/index. html 
^ http;/ /www. darkenergysurvey.org/ 
^ http://www.lsst.org/ 

http://www.euclid-ec.org/ 
^ http://wfirst.gsfc.nasa.gov/ 



Since most of the useful cosmological information con- 
tained in the cosmic shear signal is associated with small 
angular scales that are affected by nonlinear clustering, 
we also need to include these nonlinear effects to accu- 
rately model the wea k-lensing statistics (jTakada fc JainI 
12004 ISato et al.ll2009l) . Most researchers use fitting for- 
mulas based on numerical simulations or phenomenologi- 
cal approaches but it would be useful to obtain analytical 
methods that are more directly related to the cosmological 
parameters and primordial fluctuations. 

In this paper, we examine the performance of the theo- 
retical m odeling of the 3D matte r densit y distribution pro- 
posed bv lValageas fc Nishimichil ()2011al fb[). which is based 
on a combination of perturbation theories and halo mod- 
els. We focus on the convergence power spectrum and bis- 
pectrum, which are basic statistical me asurements in weak 
lensing studies (see, iHikage et al.ll20lil for a recent method 
of lensing power spectrum measurement). As compared 
with simple fitting formulas or direct numerical simula- 
tions, a significant advantage of our approach is that we 
can evaluate and compare different contributions that can 
be measured in weak-lensing surveys. Since different contri- 
butions suffer from different theoretical uncertainties, this 
is useful to estimate the accuracy that can be aimed at in 
weak-lensing statistics, as a function of scales. Furthermore, 
we find that our model provides better agreement with nu- 
merical simulations than other existent models. 

This paper is organized as follows. In Sect. [5] we first 
present our model for the 3D matter density power spec- 
trum and bispectrum. Next, we recall how this yields the 
weak lensing convergence power spectrum and bispectrum 
through the Born approximation. We describe our numer- 
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ical simulations and the data analysis in Sect. [3l Then, 
we present detailed comparisons between the simulation 
results, previous models, and our theoretical predictions, 
for the convergence power spectrum in Sect. |4l and for the 
convergence bispectrum in Sect. [Sj considering the cases of 
both equilateral and more general isosceles configurations. 
We study the relative importance of the different contri- 
butions to the power spectrum and bispectrum in Sect. [51 
arising from "1-halo" , "2-halo" , or "3-halo" terms. Finally, 
we check the robustness of our model as we vary the cos- 
mological parameters in Sect.[7]and we conclude in Sect.|Hl 

2. Analytic models 

2.1. 3D matter power spectrum and bispectrum 

Our Fourier-space normalizations for the density contrast, 
its power spectrum and bispectrum, are 



<5(x) 



dke"^''5(k). 



(^(ki)^(k2)) =,5,3(ki+k2)F(fci 



and 



(1) 



(2) 



(^(ki)^(k2)^(k3)) = SdO^I + k2 + kg) B(fcl, fc2, fcs)- (3) 

In this section we recall the model that we use to describe 
the 3D matter distribution, whic h is p resented in greater 
detail in Valageas & Nishimichi (2 011al lb[). It combines sys- 
tematic perturbation theory, which governs large scales, 
with a phenomenological halo model, which governs small 
scales. This is achieved by writing the 3D matter power 
spectrum as 



P2H+ih(A:) = P2H(A:) + PiH(fc), 



(4) 



as in the usual halo model (jCoorav fc Shethl[20nl . The 
main improvements with respect to previous works are that 
the "2-halo" term is obtained from a perturbative resum- 
mation scheme, which is complete up to one-loop order (and 
includes partial resummations of diagrams at all higher or- 
ders), while the "1-halo" term includes a "counterterm" 
associated with mass conservation that ensures its well- 
behaved asymptote at low k. 

More precisely, the 2-halo contribution is written as 



P2h(A:) = F2H(27r/fc)Pp„.t(A:), 



(5) 



where -F'2h('7) is the fraction of pairs, with initial (i.e. 
Lagrangian) separation g, that belong to two distinct ha- 
los, and Ppert{k) is the matter powe r spectrum obtained by 
perturb ation theory. As discussed in lValaeeas &: Nishimichil 
(|2011aD . it is not possible to use standard perturbation the- 
ory (beyond linear order) for Ppcrt(fc), unless one adds a 
nonperturbative ad-hoc cutoff, because this would yield a 
term that keeps growing at high k and prevents a good 
agreement with numerical simulations. A natural cure to 
this problem is to use a resumniation scheme t hat re- 
mains w ell-behaved at high k. As in I Valageas fc Nishimichil 
we u se the "di r ect st e epest- descent" resummation 
developed in IValaeeasI (|2007L l2008[) . going to "one-loop" 
order. This provides a perturbative term Ppcrt(fc) that is 
consistent with standard perturbation theory up to one- 
loop order (i.e. up to order P|) while keeping close to Pl at 



high k (thanks to the partial resummation of higher orders). 
The prefactor F2H(27r/fc) is equal to unity at all orders of 
perturbation theory, since it also writes as 



F2Yi{q) = 1-Fih((7), 



(6) 



where Pm is the fraction of pairs of initial separation q that 
belong to a single halo, and Pih = at all orders of pertur- 
bation theory (because this is a nonperturbative quantity 
that behaves as e"''^-/'^^'^^ for Gaussian initial conditions). 
Next, the 1-halo contribution is written as 



^iH(fc)^/ -m^. 



UM{kf ^W{kqM)'' 



(7) 



where f{v) defines th e halo mass function through 
(iPress fc SchechtedfTOTl 



n(M)AM = -^J{y) — , with V = 



(8) 



Here, a{qM) is the rms linear density contrast at mass M, 
or Lagrangian radius qM , with 



47r o 

and 



cr2(q)=47r/ dfc fc^P^ (fc)VF(fcq)^ 
^0 



(9) 



(10) 



where W{kq) is the Fourier transform of the top-hat of 
radius q, defined as 



Wikq)^ I ^,n..,^^ Mkq)-^kq^CO.{kq) 



(11) 



We use the scaling function /(z/) from I Valageas! (|2009t ) with 
the threshold 6l* associated with halos defined by a nonlin- 
ear density contrast of 200. In Eq.© uuik) is the Fourier 
transform of the density profile pm{x) of halos of mass M, 

UM{k) = ^—j-. = — / dxe-'^-^pA/ X , 12 

whereas W{kqM) was defined in Eg (ITTIl. We use the usual 
"NEW" halo profile (jNava rro et al.lll997D with the mass - 
concentration relation from Valage as fc Nishimichil (|2011al ) . 
The counterterm in Eq.© ensures that the 1-halo con- 
tribution decays as Pih(^) oc k^ at low fc, so that the total 
power dU goes to the linear power for CDM cosmologies. 
This follows from the conservation of matter and the fact 
that halo forma tion correspo nds to small-scale redistribu- 
tion of matter ( Peebleslll974f) . In fact, taking also momen- 
tum co nservation int o account would yield an even steeper 
k'^ tail (|Peebleslll9"7i ) but the form © is sufficient for prac- 
tical purposes. 

Thus, this framework ensures that the power spec- 
trum (|3|) is consistent with perturbation theory (up to 
the order of truncation of the resummation scheme) 
while going smoothly to the 1-halo power at high k. 
A more detailed derivatio n of E qs.(0 and ([7|) is given 
in fValagcas &: Nishimichil (|2011al) , through a Lagrangian 
point of view. 
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As described in lValaeeas &: Nishimichil (|2011b[ ). the 3D 
matter bispectrum is obtained in a similar fashion as the 
sum of 3-halo, 2-halo and 1-halo terms, 

B{ki,k2,k3) = B3ii{ki,k2,k3) + B2ii{ki,k2,k3) 

+ SiH(fci,fc2,fc3), (13) 

with 

Bm{ki,k2,k3) = -Bport(fci, ^2, fcs), (14) 



B2H(fcl,fc2,fc3) = PL{ki) J ^=7^/(^^) 



1/ p(27r) 



X n ("M(fcj) - Wik^QM)) + 2cyc., (15) 

where "2 eye." stands for two terms obtained by circular 
permutations over {fci,fc2,A;3}, and 



M 



X Y[ [uMikj) - Wik^QAi' 



(16) 



Again, the counterterms W in Eqs. dTSI) and (|16p ensure 
that the 2-halo and 1-halo contributions decay on large 
scales so that the bispectrum goes to the perturbativc pre- 



diction B, 



pert 



However, as found in Valageas & Nishimichi 



(|2011bl) and contrary to the situation encountered for the 
power spectrum, the standard one-loop perturbation the- 
ory prediction for Bpert is well-behaved at high k (i.e. it 
is significantly smaller than the 1-halo contribution) and it 
is more accurate than two alternative resummations that 
have been investigated. Therefore, we simply use the stan- 
dard one- loop perturbation theory prediction for Sport- 

While Eq. p^ yields a reasonably good match to numer- 
ical simulations (^ 10%) over all scales for the bispectrum, 
Eq.(|3]) significantly underestimates the power spectrum on 
the transition scales (by ^ 20%) even though it gives a 
good accura cy on larger scales (^ 1%) and smaller scales 
(- 10%). In IValageas fc Nishimichil (1201 Ibi l we devised a 
simple geometrical modification to Eq.Q, Ptang{k), such 
that 



-Ptang(fc) = P2H+m{k) for k <k- or k>k'+, 



(17) 



and 



log 



-Ftang(fc) 



iog(fc/fc-) 
iog(fc;/fc_) 



log 



^tang(fc4.) 



-Ptang(fc— ) 



for fc_ < k < k'j^. In other words, Ptang is iden- 
tical to P2H-1-1H of Eq.® outside of the transition 
range [fc_,fc^^], and we interpolate by a straight line in 
the plane {logfc,logP} over the interval This 
partly cures the underestimate of the power spectrum 
in this range, while keeping the perturbativc and 1-halo 
behaviors on large a n d sma ll scales. As explained in 
IValageas &: Nishimichil ()2011bf ). the wavenumbers k- and 
k'_^ are determined by the constraint that the reduced equi- 
lateral bispectrum, Qcq{k) = B{k,k, k)/[3P{k)'^], is mono- 
tonically increasing. Indeed, it happens that Eq.Q under- 
estimates the actual power spectrum on transition scales, 



which leads to a spurious peak for Qcq{k). Then, from the 
shape of this naive prediction for Qcq{k) we automatically 
detect the range , k'^] where the model is not sufficienly 
accurate (because higher-order perturbative contributions 
come into play or the decomposition over 2-halo and 1- 
halo terms is too simple to accurately describe regions 
such as outer shells and filamentary branches of collaps- 
ing halos that have not yet relaxed). This typically gives 
fc_ = 0.4/iMpc"i and fc+ = 3/iMpc"^ at 2; = 1. The main 
point here is that the range , k'_^_] is automatically de- 
termined from the shape of P2H+mik) and B{k, fc, fc), and 
it evolves with r edshift and with the in i tial po wer. We re- 
fer the reader to IValageas fc Nishimichil (j2011bt ) for further 
details. 

Thus, Eqs.(IT3|)-(ITn|) and 1^7^-^ define the analytic 
model that we investigate in this paper for the 3D matter 
bispectrum and power spectrum. 

2.2. 2D weak-lensing power spectrum and bispectrum 

Using Born's approximation, the weak-lensing convergence 
k{6) can be written as th e integral of the density con- 
trast along the line of sight (jBartelmann fc Schneideiir200ll : 
IMunshi et aIll2008D . 



Jo 



(19) 



where x and 2? are the radial and angular comoving dis- 
tances. 



w{x,Xs) 



3n^,H^v{x)v{xs - x) 

2c^V{xs) 



(1 + z), 



(20) 



and Zs is the redshift of the source (in this article we only 
consider the case where all sources are located at a single 
redshift, to simplify the comparisons with numerical simu- 
lations and the dependence on the source redshift). Then, 
using a flat-sky approximation, wh ich is valid fo r small an- 
gles below a few degrees ( Valageas et al.|[2"01lD . we define 
its 2D Fourier transform through 



<e) = J 



= I d£e'^-^^(£). 



(21) 



As in 3D, we define the 2D convergence power spectrum 
and bispectrum as 

{~K{£i)k{£2)) =SD{il+£2)Pn{il), (22) 

and 

{ii{ei)k{e2)Kie3)) = 5d{£i + £2 + £3) B^{£i,£2,i3). (23) 

(18) From Eq.^ on e obtain s at on ce, using Lir nber's 



approximation (iLimbeil 119531: iKaiseil 119921: 



iBartelmann fc Schneideill200ll : iM^lnshi et all 1 20081) . 

PM^^^ r^X^Pie/V-^z), (24) 



B^{£i,£2,i3) = (27r) 



2 1^ dx^Bi£,/V,£2/V,£3/V;z).i25) 



This provides the weak-lensing convergence power spec- 
trum and bispectrum from our model described in Sect. 12. ll 
for the 3D matter density power spectrum and bispectrum, 
through a simple integration over the radial coordinate up 
to the source plane. 
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3. Numerical simulations 

To obtain accurate predictions of the convergence power 
spectrum and bispectrum we must perform the ray-tracing 
simulations through large-volume, hig h-resolution TV- 
body simulations of stru c ture formation (I Jain et alJ 20001: 
Hamana fc Mellied l200lt iHilbert et all l2009t ISato et alJ 
2009() . To perform the TV-body sirnulations, we use a modi- 
fied version of the Gadget-2 code (jSprin gcl 2005). The ray- 
tracing simulations are constructed from 2 x 200 realizations 
of A^-body simulations with cubic 240 and 480/i~^Mpc on 
a side, respectively. We employ 256"^ particles for each sim- 
ulation. For our fiducial cosmology, we adopt the standard 
ACDM cosmology with matter fraction f^m = 0.238, baryon 
fraction fib = 0.0416, dark energy fraction J7de = 0.762 
with the equation of state parameters wq = —I and Wa — 0, 
spectral index Ug = 0.958, normalization As = 2.35 x 10~^, 
and Hubble parameter h = 0.732, which a.re con sistent with 
the WMAP 3-year results (jSoergel et al.ll2007D . This fidu- 
cial cosmology gives for the variance of the density fluc- 
tuations in a sphere of radius 8/i~^Mpc the normalization 
(78 = 0.759. We consider source redshifts at either Zg — 0.6, 
1.0, or 1.5. Using ray-tracing simulations we generate 1000 
realizations of 5° x 5° convergence maps for each source 
redshift. 

In addition to the fiducial cosmology case, we performed 
ray-tracing simulations for several slightly different cos- 
mologies. We vary each of the following cosmological pa- 
rameters: As, Us, the cold dark matter density Vl^h? , ^dc, 
and Wq by ±10%, respectively. In varying the parameters, 
we keep the flatness of the universe and the physical baryon 
density ^\^h? to be unchanged. Therefore the three param- 
eters, h, f2m, and f^b, are varied simultaneously. For each 
of these 10 different cosmologies, we obtain 40 realizations 
of convergence fields for each of the three source redshifts. 
Details of the me thods used for the ray-tracing simulations 
can be found in ISato et all (|2009D . All converg ence maps 
used i n this paper are the same as those used in lSeo et al.l 
(|201lh . but for varied cosmologies we use three orthogonal 
projection axes instead of using only one projection axis to 
increase the independence. 

In Sects. 010 we use the maps for the fiducial cosmology 
while in Sect. [7] we use those for the varied cosmologies to 
investigate the robustness of our model. In Sect. [7] we show 
the results for six cases, varying rt^, ^Ich^ , and wq by ±10%, 
which we denote as ns±, Vlch?±, and u'o±. The exact values 
for these cosmological parameters are listed in Table [XT] in 
App.El 

The binned convergence power spectrum and bispec- 
trum are measured as follows. We first apply Fast Fourier 
Transformation to each of the convergence fields to obtain 
k{l). We then bin the data into logarithmically-equal bins 
in ^, whose widths are set as Aln£ = 0.3 for the power 
spectrum, and Aln^ = In 2/2 « 0.35 for the bispectrum, 
respectively. The power spectrum and the bispectrum are 
obtained by simply averaging the products of modes: 



(26) 



B«(4,4,4) = V Re[k{l^)k{l2)h{l^)], (27) 

where Re[...] denotes the real part of a complex number, and 
the summation runs over modes t. (jEi, i — 1,2, 3) which fall 



into bin i (ii) for the power (bi-) spectrum. In the above, 
N( is the number of modes taken for the summation and 
is given by Ng « ^shcui^s/(27r)^ where Aghcn is the area of 
two-dimensional shell around the bin £ and can be given 
as ^shoii ~ 2tt£A£ + tt{A£)'^, and ils is the survey area. 
Similarly, the factor iV^^ , which appears in the estima- 
tor of the bispectrum, denotes the number of triangles in 
£ space. We then average the measured spectra over 1000 
random realizations to obtain our final estimates of Pk and 
B^. The statistical uncertainty of our estimates are also 
estimated from the variance of and over 1000 real- 
izations. We plot the 3-cr uncertainty as error bars in what 
follows. 

4. Convergence power spectrum 

We now compare our results for the weak-lensing power 
spectrum with numerical simulations. In addition, we also 
consider the predictions obtained from the popular "halo- 
fit " fitting funct i on for the 3D density power spectrum given 
in ISmith et al.l ()2003D . to estimate the gains that can be 
reached by using a somewhat more systematic approach. 

We show our results for the convergence power spec- 
trum, Pk{£), in Fig. [1] For reference we also plot the linear 
power and we can check that both simulations and analyt- 
ical models recover the linear regime on large scales. The 
numerical error bars increase on large scales because of the 
finite size of the simulation box. On small scales the nu- 
merical error is dominated by systematic effects, due to 
the finite resolution, that lead to an underestimate of the 
small-scale power (as clearly shown by the sharp decline at 
very high £) . For each source redshift we estimated the scale 
£ up to which the simulations have an accuracy of better 
than 5% by comparing with higher-resolution simulations 
(with 512'^ particles instead of 256"^). This scale is shown 
by the vertical arrow in Fig. [T] and we can check that our 
model indeed agrees with the numerical simulations up to 
this multipole. 

We recover the well-known fact that on scales that 
are of interest for cosmological analysis of weak gravi- 
tational lensing, 10^ < £ < lO'*, nonlinear terms sig- 
nificantly contribute to or dominate the pow er spectrum 
(|Bartelmann fc SchneideJ 120011 : llvkinshi et al.l [20081. The 
overall shapes of the power spectra obtained from our 
mod el, described in Sec t. [H and from the fitting formula 
from lSmith et al.l (|2003l ). are very similar. This is expected 
since the latter "halo-fit" also goes to the linear power on 
large scales and it is based on numerical simulations for its 
small-scale behavior, while the halo model that underlies 
our approach on small scales also agrees with similar sim- 
ulations (in particular it involves the "NFW" profile and 
a mass-concentration relation that are derived from simu- 
lations). These similar shapes also confirm that the sharp 
falloff of the power at high £ found in the ray-tracing sim- 
ulations is not physical but due to the finite resolution. 

The improvements that can be expected from our ap- 
proach, as compared with a simpler fitting formula to mea- 
sures of the 3D power, are that i) on large scales we are 
consistent with perturbation theory up to one-loop order, 
and ii) on small scales we directly use a physical halo model 
(instead of using some formula for the power spectrum with 
free exponents that are fitted to a set of simulations). 

On large scales the error bars of our ray-tracing simula- 
tions (and their lack of large-scale power) are too large to 
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Fig. 1. Convergence power spectrum for sources at redshifts Zg = 0.6, 1, and 1.5. The points are the results from numerical 
simulations with 3 — cr error bars. The low black dash ed line "L" is the linear power, the middle blue dot-dashed line "S" 
is the result from the "halo-fit" of ISmith et al.l ()2003l ). and the upper red solid line "comb." is the result from our model, 
which combines 1-loop perturbation theory with a halo model. The vertical arrow close to the peak shows the scale up 
to which the simulation result is valid within 5%. 



clearly see the benefit of the one-loop perturbative terms at 
low Zg. However, the higher accuracy due to these higher- 
order perturbative contributions can be seen at = 1 and 
Zs = 1.5, as will be shown more clearly in Sect. 16.11 be- 
low. (The benefit of these terms has already been shown 
in studies of the 3D matter density power spectrum, espe- 
cially for the a ccurate prediction of the baryon acoustic os- 
cillati o ns, e.p . iJeong fc Komatsul (l2006|): Nishimichi et al.l 

2008"); 



20071 mo. 



2008) 



Crocc e fc Scoccimarro' 

Taruva et al.. (20091. Sato & Matsubara (2011); 

IValaeeas fc Nishimichi! (|2011al lb[)). In the weak-lensing con- 
text, this higher accuracy could also be useful for the 
analysis of future obs ervations such as the Euclid mission 
(|Refregier et al.ll2010D . 

On small scales, although the agreement of our model 
with the simulations is not perfect we can see a clear 
improvement as co mpared with the fitting formula from 
[Smith et al^ ([200^. This is not surprising since the lat- 
ter was derived from a set of older numerical simulations 
with somewhat different cosmological parameters than the 
ones we consider in this paper. The underestimate of the 
power on small scales by this "halo-fit" formula, by about 
30% around t he pe ak for Zg = 1, was already noticed in 
iHilbert et all (120090 . The comparison of our model wi th 
the 3D power spectrum in I Valageas fc Nishimichil ()2011al lbl) 
was also based on simulations with different cosmological 
parameters, but through the explicit expression ([7]) we au- 
tomatically take into account the dependence on cosmol- 
ogy of the linear growth factor D^{t) of the density fluc- 
tuations and of the linear threshold 5l* associated with 
the nonlinear density contrast of 200 (Valageas 2009). We 
mostly neglect the dependence on cosmology of the halo 
profiles and of their mass-concentration relation. However, 
even though virialization processes do not reach complete 
relaxation, that is, a violent relaxation "a la Lynden-Bell" 
(jLvnden-Belll 119671 ) is not complete in this cosmological 
context and does not allow the halos to reach a statisti- 
cal equilibrium that is fully independent of the initial con- 
ditions, internal halo properties are almost universal and 
independent of cosmological parameters up to a good accu- 
racy (for realistic CDM scenarios). Similarly, deviations of 
the halo mass function from "universality" have been de- 



tected but_jx£ra;ther_wrea iMore et all 

1201 ll : iBhattacharva et al.ll2011 ). This is especially true for 
our purpose as we integrate over the full halo mass function 
(which is normalized to unity) and we take into account the 
dependence on cosmology of its large-mass cutoff. This ex- 
plains why our model is able to reach a good agreement 
with the numerical simulations. 

On the intermediate scales, which are sensitive to the in- 
terpolation ([TS]) , we also obtain a satisfactory match to the 
numerical simulations. Within this range, around i ^ 300, 
the "halo-fit" prediction happens to provide a similar, or in 
a few cases slightly better, agreement with the simulations. 
These transition scales, where both the one-loop pertur- 
bative contribution and the 1-halo term are subdominant 
as shown in Fig. 2] below, are governed by the interpola- 
tion (|18p. This means that this interpolation is not fully 
satisfactory and that there remains room for improvement. 
However, to keep our model as simple as possible and to 
remain consistent with our previous 3D studies we do not 
investigate here alternative prescriptions. 

We shall check in Sect. [7] below that the overall agree- 
ment found in Fig. [T] remains valid as we change the values 
of the cosmological parameters. 



5. Convergence bispectrum 

We now compare our results for the weak-lensing bispec- 
trum with numerical simulations. In addition, we also con- 
sider the predictions obtained from the following four sim- 
ple models, which have been used in some previous works. 

We first compute the results obtained from the lowest- 
order ( "tree-order" ) prediction from sta ndard perturba tion 
theory for the 3D bispectrum (Bernard eau et aDl2002l) . 

Bt,,,^{ki,k2,h)^2F2{ki,k2,fii2)PL{ki)PL{k2)+2cyc.{28) 

where /ii2 = (ki • k2)/(fcifc2) and 



^2(^1,^2,^^12) = 



1 / fci k2 



7 2 \k2 ki 



A*12 



7 



Ml2- 



(29) 
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Fig. 2. Convergence bispectrum for sources at redshifts Zg = 0.6, 1, and 1.5, for equilateral triangles. The points are the 
results from numerical simulations with 3 — cr error bars. The low black dashed line "trecL" is the tree-order bispectrum 
(P5|) . t he green dot t ed lin e "treeNL" is the tree-order formula (PH)) . where we use the nonlinear power from the "halo- 
fit" of ISmith et al.' (2003"), the two yellow dashed lines "F2.nl" (I51T) include in addition a fitting formula for the kernel 
-p2,NL from [Scoccimarro & Couchman (2001) and use either Ps{k) (lower curve) or Ptang(fc) (upper curve) for the 3D 
power spectrum, the upper blue dot-dashed line "seal." is the phcnomenological scale transformation ([5^ from lPan et al.l 
(j2007() . and the red solid line "comb." is our combined model (fT5|) - ([T5| . The vertical arrows are at the same scale as in 
Fig.ffl 



Second, we consider a "tree-nonlinear" approximation 
where in Eg. ((28)) we replace the linear 3D power Pj ,(k) by 
the nonlinear power Ps{k) from ISmith etaP (l2003l) . 

Btroo«L iki,k2,k3)^ 2F2{ki,k2, ^il2) Psiki)Ps{k2)+2 eye. (30) 

Third, we consider the fitting formula from 
IScoccimarro fc CouchmanI ()200lD . 

BF2,NL{ki,k2,k3) = 2F2,-Mhiki,k2,Hi2)Ps{ki)Ps{k2) 

-f2cyc. (31) 

where the kernel -F2,nl is an effective kernel that inter- 
polates from the large-scale perturbative result (1^ to a 
small-scale ansatz where the angular dependence on /ii2 
vanishes. He re we shall u s e the nonlinear power spectrum 
Ps{k) from ISmith et al.l (|2003ll as well as the nonlinear 
power spectrum Ptang(fc) of our model, Eas. (fT7| - (fT8)) . 

Four, we con sider the scale transformation studied in 
iPan et aP (|2007l) . foll owing the spirit of the s cale transfor- 
mation introduced in iHamilton et al.l (Il99lt) for the two- 
point correlation function and next in iPeacock fc Doddsl 
(|1996[) for the power spectrum, 

Bscal(fcl,fc2,fc3) = w{ki)w{k2)w{k3)2F2{ki,k2,fj-i2) 

xPL(fci)PL(fc2)+2cyc. (32) 

with for i = 1, 2, 3, 

k^[i + Al{h)]-^/^h, wCh) - ^Ps{h)/PL{k). (33) 
5.1. Equilateral triangles 

We first show in Fig. [2] our results for the bispectrum for 
equilateral configurations. By construction, all curves must 
converge to the lowest-order perturbative prediction (^51) on 
large scales. We can check that this is indeed the case. As for 
the lensing power spectrum our simulations are too small 
to reach this "tree-order" regime but they are consistent 
with this asymptotic behavior. 



Again, we recover the well-known fact that on the scales 
of interest lowest-order perturbation theory is not suffi- 
cient to predict weak-lensing statistics and we must take 
into account higher-order terms or nonperturbative con- 
tributions. Using the nonlinear power within the "tree- 
level" expression, as in Eg. ([50)1 . significantly increases the 
convergence bispectrum and improves the overall shape 
but it remains far from the simulation results and can- 
not be used for practical purposes. Adding a modification 
to th e kernel F2, as in Eg. dSTl) using the fitting formula 
from IScoccimarro fc CouchmanI ()200l[ ). further improves 
the predictions and it provides a good match to the sim- 
ulations on weakly nonlinear scales. However, this simple 
fitting procedure underestimates the power at high t. This 
is partly related to the underestimate of the power spec- 
trum on small sca les by the "halo-fit" fitting formula of 
ISmith et al.l (|2003D . Indeed, we can see that using in Eq. pip 
the nonlinear 3D power spectrum provided by our model, 
which provides a good match to simulations for both the 
3D power ([Valageas &: Ni shimichi 2011b) and the 2D con- 
vergence power spectrum (Fig. |T]) increases the power at 
high I. This improves the agreement with the simulations 
but there remains a significant discrepancy. 

The scale transformation (j32l) fares better on weakly 
nonlinear scales (at least for this equilateral configuration) 
but it greatly overestimates the power at high I. The same 
breakdown was already observed for the 3D bispectrum 
(|Pan et al.ll2007l: I Valageas fc Nishimichill2011b[ ). The best 
agreement with the numerical simulations is provided by 
our model (fT5|) - (IT5|) . In particular, it is interesting to note 
the good match on the transition scales, i ~ 10^, which a 
priori are the most difficult to reproduce since they are at 
the limit of validity of both perturbative approaches (which 
break down at shell crossing) and halo models (which as- 
sume virialized halos) . A similarly g ood agreement was also 
observed for the 3D bispectrum in IValageas fc Nishimichil 
(l20TTbh . 

At high i we again predict more power than is measured 
in the simulations, but as for the lensing power spectrum 
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Fig. 3. Convergence bispectrum ^^(^i, ^2, ^3) for isosceles configurations, where £1 — £2- The first two rows show the 
dependence on the common side £1 = £2, at fixed £3 (in the weakly nonlinear regime in the first row, and in the highly 
nonlinear regime in the second row), while the third row shows the dependence on the third side £3 at fixed £1 = £2 
(around the transition scales). The symbols are the same as in Fig. [5] 



this is at least partly due to the lack of small-scale power 
in the simulations because of the finite resolution. Thus, we 
again plot in Fig. [2] the vertical arrows that were plotted in 
Fig.[T] Since the bispectrum typically scales as the square of 
the power spectrum this should roughly correspond to an 
accuracy threshold of about 10% for the simulations. We 
can check that our model agrees with the numerical results 
up to this multipole. 

5.2. Isosceles triangles 

We now plot in Fig.|3]the convergence bispectrum for isosce- 
les configurations, where £1 = £2- Instead of the angle ■di2 
between the two sides £1 = £2 we show the dependence on 



the length of either the two equal sides or of the third side. 
This avoids putting all "squeezed" triangles in a narrow 
region around ??i2 = 0. 

As for the equilateral configurations shown in Fig. [2J we 
can check that the "tree-order" bispectrum is not the 
dominant contribution on these scales, inserting the non- 
linear power spectrum as in (1301) is not sufficient, and the 
scale transformation p2p breaks down too early on nonlin- 
ear scales. Thus, the only two reasonable models are the 
expression pil) . which involves fitting formulas for both 
the nonlinear power spectrum and the kernel i^2,NL, and 
our model ((i3l)-(Il6l). 

In the first row, where £3 ~ 107 (the analysis uses 
the fiat-sky approximation hence £ is not necessarily an 
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integer), both curves (j31l) and (IT3l) agree rather weU 
with the numerical simulations, except that for extremely 
"squeezed" configurations, £i ^ lOO^a, our model yields 
a downturn which underestimates the bispectrum whereas 
the fitting model pip remains in good agreement with 
the numerical data. We shall come back to this point in 
Sect. 16. 2T2I below, where we shall find out that this is prob- 
ably due to an inaccuracy of the 2-halo contribution -B2H- 
On the other hand, using the nonlinear power spectrum 
-Ptang(^) in Ea. ((3T|) increases the power on small scales, as 
in Fig. [5J which spoils the match to simulations for these 
extremely "squeezed" configurations. 

In the second row, which is farther into the nonlinear 
regime, we find in agreement with the equilateral case of 
Fig. [2] that the fitting formula ([5T|) underestimates the bis- 
pectrum whereas our model provides a reasonable match to 
the numerical simulations. More importantly, Ea. (l3ip does 
not reproduce the shape of the bispectrum as a function 
of £1, since it does not capture the falloff that appears for 
£1 ^ 10^3, and this is not cured by using a more accurate 
3D power spectrum. In contrast, our model (jl3p correctly 
predicts this shape as well as the overall amplitude. 

In the third row, which corresponds to transition scales 
in the mildly nonlinear regime, neither model shows a per- 
fect match to the simulations but our prediction fares sig- 
nificantly better and still provides a reasonable agreement. 

These results show that it is difficult to reproduce the 
bispectrum for a wide variety of configurations and scales 
by using a global ansatz such as (|3T|) . Approaches such as 
the one studied in this paper, which are built on explicit 
and physical ingredients, are better controlled and offer an 
easier route to systematic improvement. Indeed, by improv- 
ing the accuracy of each ingredient (e.g., including higher 
orders of perturbation theory or more complex halo pro- 
files) one should reach a higher accuracy for the final 3D 
and 2D statistics for any configuration and scale (although 
transition regimes may still prove difficult). Moreover, by 
splitting the problem into several elements it is easier to 
make progress by improving each contribution in turn. 

6. Relative importance of the different 
contributions 

We have seen in the previous section that our model, based 
on a combination of perturbation theories and halo mod- 
els, provides a good match to numerical simulations. This 
means that such an approach can be useful to obtain predic- 
tions for weak-lensing statistics for a variety of cosmologies, 
which is an important goal for observational and practical 
purposes. 

A second advantage of our approach, as compared with 
simple fitting formulas (or direct simulations), is that we 
can easily evaluate and compare the different contributions 
that eventually add up to the signal that can be measured in 
weak-lensing surveys. Thus, we can distinguish the various 
perturbative terms as well as the nonperturbative contri- 
butions associated with 1-halo or 2-halo terms. 

This is useful to estimate the accuracy that can be 
aimed at in weak-lensing statistics, as a function of scales, 
since different contributions suffer from different theoretical 
uncertainties. For instance, the perturbative contributions 
are obtained in a systematic and rigorous manner for any 
cosmology, and the accuracy is only limited by the order of 



the perturbative expansion, whereas "halo contributions" , 
which depend on the halo profiles and mass functions, are 
phenomenological ingredients that involve at some stage in- 
put from numerical simulations and do not offer systematic 
routes to arbitrarily high accuracy. 

This also allows us to understand which aspects of the 
dynamics are probed by weak-lensing statistics, so that de- 
pending on one's goals (e.g., to constrain cosmological pa- 
rameters or to estimate halo properties themselves) one can 
focus on the appropriate range of the weak-lensing signals. 

6.1. Convergence power spectrum 

We plot our results for the convergence power spectrum in 
Fig. m showing the underlying 2-halo and 1-halo contribu- 
tions in addition to the full model curve that was already 
shown in Fig. [T] We can see that including the 1-loop per- 
turbative term yields slightly more power on quasi-linear 
scales and helps to obtain a good match to the simula- 
tions. At high £ the 2-halo contribution goes back close 
to the linear power thanks to the partial resummation of 
higher orders. As recalled in Sect. 12.11 this is a useful im- 
provement over the standard 1-loop perturbation theory 
because it ensures that at high k for the 3D power, and 
at high £ for the lensing power, the 2-halo term does not 
give significant contributions. At high £ around its peak, 
the lensing power becomes dominated by the 1-halo term 
but as for the 3D power there remains a wide intermediate 
range, about a decade over £. In our model, the 3D power on 
these transition sca les is described by the iri t erpolat ion (fT8|) . 
which was seen in IValageas fc Nishimichil ()2011bD to pro- 
vide a good match to numerical simulations. As expected. 
Fig. |4] shows that this also provides a good interpolation 
for the 2D convergence power spectrum. It may be possible 
to reduce this range on its low-^ side by including higher 
orders of perturbation theory in the 2-halo term, but it is 
likely that the high-^ side receives significant contributions 
from nonpertur bative terms, associated with shell crossing 
([Valageas! l201l[ ) , that are not well described either by a 
simple 1-halo term (e.g., associated with dense filaments 
or infalling regions close to the virial radius). In any case. 
Fig. |4] clearly shows how the convergence power spectrum 
depends on large-scale perturbative density fluctuations or 
on small scale halo properties, as the multipole £ varies. 

Future surveys plan to measure the weak lensing power 
spectrum on these scal es, 50 < £ < 7000, with an accuracy 
of about 1% for LSST (ILSST Science Collab oration s et al. 
200^ or even slightly better for Euclid (Laureii s et al. 
12011 1. Our numerical simulations do not allow us to check 
our model down to this accuracy. However, the comparison 
between our predictions and the "halo- fit" in Fig. [T] already 
shows that this is beyond the accuracy of current models 
and is a very challenging theoretical goal. Our model al- 
ready provides an accuracy of about 1% on perturbative 
scales and 10% on the transition and highly nonlinear scales 
(this is the typical acc u racy ch ecked on 3D simulations in 
IValageas fc Nishimichil ()2011bD ). Figure 01 explicitly shows 
that for 100 < £ < 1000 the one-loop correction to the lin- 
ear power is required to reach this accuracy. Unfortunately, 
for £ > 200 this is not sufficient since higher-order terms 
or non-perturbative corrections cannot be neglected and 
the power spectrum is sensitive to the transition scales, de- 
scribed by the interpolation (IT51) . This suggests that for 
£ > 500 a percent accuracy is probably beyond the reach of 
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Fig. 4. Convergence power spectrum for sources at redshifts Zs = 0.6, 1, and 1.5. The points are the results from numerical 
simulations with 3 — a error bars. The low black dashed line "L" is the linear power, the middle black dot-dashed line 
"1-loop" is the result from the 2-halo contribution ([S]), where we use a perturbative resummation that is complete up to 
1-loop order, the upper blue dashed line "IH" is the result from the 1-halo contribution ([S]), and the red solid line is the 
resuh of our model (fT7 | - (fT8)) . 




Fig. 5. Convergence bispectrum for sources at redshifts Zs — 0.6, 1, and 1.5, for equilateral triangles. The points are the 
results from numerical simulations with 3 — a error bars. The low black dashed line "L" is the tree-order bispectrum 
([25]) . the "1-loop" dot-dashed black line is the prediction of 1-loop standard perturbation theory, which is identified with 
our 3-halo term p^ . the green dotted line "2H" is the 2-halo contribution (IT5|) and the upper blue dashed line "IH" is 
the 1-halo contribution ([TB| . 



systematic perturb ative approache s, which cannot go be- 
yond shell crossing (' Valageaj|201lD , and requires dedicated 
numerical simulations. 

6.2. Convergence bispectrum 

We now study the contributions from various terms to the 
lensing bispectrum. From Ea. (ll3p we now have three con- 
tributions, associated with the 3-halo, 2-halo and 1-halo 
terms, the 3-halo contribution being identified with the per- 
turbative contribution. 

6.2.1. Equilateral triangles 

We plot our results for the convergence bispectrum for equi- 
lateral configurations in Fig. [5l showing the various con- 
tributions in addition to the full model curve that was 
already shown in Fig. [5] As recalled in Sect. 12.11 the 3- 
halo term is identified with the perturbative prediction and 
in this paper we use standard perturbation theory at one- 



loop order for this contribution . As for the 3D bispectrum 
(jValageas fc Nishimichi|[2011bl) . and contrary to the power 
spectrum, this gives a contribution that becomes negligi- 
ble on small scales, so that it is not necessary to use a 
resummation scheme or to add a nonperturbative cutoff to 
ensure a good high-^ behavior. On the other hand, con- 
trary to the power spectrum we can see that going to 1- 
loop order now makes a strong improvement over the tree- 
order result. This feature was already noticed for the 3D 
density field (jSefusatti et al.l[2010t IValageas fc Nishimichil 
2011b). In fact, combining this 1-loop perturbative con- 
tribution with the 2-halo and 1-halo terms is sufficient to 
obtain a good match to the simulations, from the quasi- 
linear to the highly nonlinear scales. This suggests that 
higher orders of perturbation theory do not significantly 
contribute to the bispectrum and that we already have a 
reasonably successful model. For these equilateral triangles 
the 2-halo term is also subdominant on all scales (by a fac- 
tor ~ 10 at least). This is a nice property since being a 
mixed term, which involves both large-scale halo correla- 
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Fig. 6. Convergence bispectrum -Bk(£i, £2, -^3) for isosceles configurations, as in Fig. [31 The symbols are the same as in 
Fig. [SJ In the first row we show the absolute value of the 1-halo contribution, which becomes negative at £1 ~ 2500. 



tions and internal halo structures, it may be more difficult 
to predict than the 3-halo term (derived from systematic 
perturbation theories) and the 1-halo term (that only de- 
pends on internal halo profiles and mass function). These 
vario us features were also observed for the 3D bispectrum 
(jVala geas & Nishimichi 2011b), from which they directly 
derive through the relation ((25)) . 

Again, our numerical simulations do not allow us to 
reach the pe rcent accur acy aimed at by future surveys 
(ILSST Scienc e Collaborat ions et al.l l2009t iLaureiis et al.l 
|2011|) . and Figs. [5] and [5] show that this is a challeng- 
ing goal. As noticed above, including 1-loop corrections is 
clearly important to reach this accuracy. A nice feature is 
that without introducing any interpolation on transition 
scales we already have a reasonably good model and that 
the 2-halo term contributes to about 1% or less on most 



scales. Then, our model already provides an accuracy of 
10 % (this is the typical ac c uracy c hecked on 3D simulations 
in IValageas fc Nishimichil ()2011b[ )) on the scales shown in 
Fig. [5J and better at low £ in the perturbative regime. At 
high £, £ > 10^, including the effects of baryons on halo 
profiles (as compared with dark matter only N-body sim- 
ulations) is probably necessary to reach percent accuracy, 
but this is beyond the scope of this paper. 

6.2.2. Isosceles triangles 

We turn to isosceles configurations in Fig. [SI As for the 
equilateral case plotted in Fig. [5l the upper row shows that 
on weakly nonlinear scales the 1-loop contribution makes a 
significant improvement over the tree-order prediction and 
allows us to obtain a good match to the numerical sim- 
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ulations. Moreover, we can see that our underestimate of 
the bispectrum for very "squeezed" triangles, at £i ^ 100^3 
with £3 ~ 107, which was aheady noticed in Fig. |31 corre- 
sponds to a change of regime. Indeed, this downturn fohows 
the 1-loop prediction and it occurs at a point where the 
2-halo and 1-halo contributions become important. Then, 
this discrepancy may be understood as a difficuhy of our 
simple model to reproduce the bispectrum in this corner of 
configuration space because the transition from the 3-halo 
perturbative regime to the 1-halo highly nonlinear regime 
is not sufficiently well described (e.g., simple spherical halo 
models cannot be expected to describe very well interme- 
diate regions such as filaments and outer infalling regions 
of halos). Alternatively, it may happen that adding higher 
orders of perturbation theory to the 3-halo contribution 
improves the agreement with simulations for these isosceles 
triangles, without contributing much to the equilateral tri- 
angles where the agreement is already rather good. Another 
possibility is that our 2-halo term is not very accurate and 
gives too small a contribution for these isosceles configu- 
rations. This would be a natural and simple cure, because 
as seen in Fig. [5] this 2-halo term is negligible for equilat- 
eral cases. Therefore, it should be possible to "tune" the 
2-halo term so as to match the isosceles simulation results 
of Fig. |6] without damaging the good agreement obtained 
for the equilateral case in Fig. [S) This would also be a more 
natural explanation than a problem in the 1-halo term, be- 
cause the latter should be simpler to model (since by defini- 
tion it only deals with inner regions of halos) and it provides 
successful results in regimes where it dominates the bispec- 
trum, as shown by the second row in Fig. [6] and the high-£ 
part of Fig. [S] 

Since the problem found in the high-^i part of the up- 
per row of Fig. ini corresponds to very "squeezed" triangles, 
£1 ^ 100^3, hence to a limited region of the configuration 
space for the triplet {£i,£2.,£3}, and our goal is to study 
simple and systematic models (i.e., that do not require fit- 
ting internal parameters for each new set of cosmological 
parameters), we do not investigate in this paper how to 
improve the model (e.g., through added complexity of the 
2-halo term) to better describe this regime. 

We can see in the second row, which is fully dominated 
by the 1-halo contribution, that the curved shape of the 
bispectrum does not arise from a change of regime and a 
switch from one contribution to the other as in the first 
row. This is a genuine effect associated with the behavior 
of the 1-halo term that is well described by our model. In 
contrast, as seen in Fig. [31 this shape was not well recovered 
by the ansatz (|3ip. This shows the benefits of our approach 
that explicitly combines perturbation theory with a halo 
model, taking advantage of the good behavior of these two 
ingredients in the regimes where they dominate. 

In the third row, which corresponds to transition scales, 
we can see a change of regime as the bispectrum is dom- 
inated by the 1-loop perturbative contribution at low £3 
and by the 1-halo contribution at high £3. In contrast to 
the first row, the 2-halo contribution does not play a signif- 
icant role and our model provides a good agreement with 
the simulations on all scales. This is a nice result since one 
could have expected transitions to be difficult to reproduce. 
This again suggests that the mismatch found for extremely 
"squeezed" triangles in the first row may be due to the 2- 
halo term. Moreover, this shows that a good modeling of 
the 3-halo (i.e. perturbative) and 1-halo terms is sufficient 
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Fig. 7. Convergence power spectrum for sources at redshift 
Zs = 1, for the six cosmologies given in Table lA.ll The 
points are the results from numerical simulations and the 
lines are the predictions of our model. 
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Fig. 8. Convergence bispectrum for sources at redshift Zs — 
1, for the six cosmologies given in Table ET] The points are 
the results from numerical simulations and the lines are the 
predictions of our model. 

to provide accurate predictions for the lensing convergence 
bispectrum for a wide variety of configurations including 
some transition regimes. This agrees with the good agree- 
ment obtained for equilateral configurations, which are also 
dominated by the 3-halo and 1-halo terms as seen in Fig.[5l 

7. Dependence on cosmology 

In this section we investigate the robustness of our model 
as we vary the cosmological parameters. We consider six 
alternative cosmologies, where n^, ^ch^, and wg are mod- 
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ified by ±10% with respect to the fiducial cosmology used 
in the previous sections. The values of the associated cos- 
mological parameters are given in Table ETTl in App. 1X1 We 
compare the predictions of our model with numerical sim- 
ulations for these six alternative cosmologies in Figs. [7] and 
[51 for the convergence power spectrum and bispectrum at 
Zg — 1. To avoid overcrowding the figures we do not plot 
the error bars of the numerical simulations. Each pair n^, 
rich^, and wq gives two curves that are roughly symmet- 
ric around the fiducial cosmology result, since we consider 
small deviations of ±10%. The deviations are largest for 
the Us case, which changes the shape of the initial power 
spectrum as well as the normalization erg. These six cases 
roughly cover the range that is allowed by current data and 
the Hs c ase is already somewh at beyond the observational 
bounds (jKomatsu et al.l[201lh . Thus, they provide a good 
check of the robustness of our model for realistic scenarios. 

We can see that the dependence of the convergence 
power spectrum on the cosmological parameters is well re- 
covered by our model. The agreement is somewhat less ac- 
curate for the bispectrum, since as for the fiducial cosmol- 
ogy the theoretical predictions give more power on small 
scales than is measured in the simulations. However, we re- 
cover the general trend and obtain a reasonable match up 
to intermediate scales. We can see that the simulation re- 
sults for the convergence bispectrum are scattered around 
the theoretical predictions on large scales because the num- 
ber of realizations is not large enough to converge. (Here we 
only use 40 realizations, whereas for the fiducial cosmology 
we used 1000 realizations.) On small scales the effects of the 
finite resolution are larger than for the power spectrum and 
the range of wavenumbers where simulation results are re- 
liable is clearly narrower in Fig. [51 than in Fig. [71 (we clearly 
see unphysical deviations at low and high i) . This explains 
why the agreement appears worse for the bispectrum as this 
is partly due to the larger numerical noise. 

We obtain similar results for Zg — 0.6 and Zs — 1.5, 
as well as for other cosmologies where we vary As or fido 
by ±10%. This shows that our model and, more generally, 
models based on combinations of perturbation theory and 
halo models provide a good modeling of the matter distri- 
bution and of weak gravitational lensing effects and capture 
their dependence on cosmology. In particular, Figs. [71 and 
ISl show that the accuracy of our model is sufficient to con- 
strain Us, flch'^, and wq to better than 10%. 

8. Conclusion 

In this article we have investigated the performance of 
current theoretical modeling of the 3D matter density 
distribution with respect to weak-lensing statistics, fo- 
cusing on Fourier space statistics, specifically the con- 
vergence power spectrum and bispectrum. We find that 
for both quantities a niodel in troduced in previous works 
([Valageas fc Nishimichill20"llal lbh that combines 1-loop per- 
turbation theory with a halo model fares better than some 
other recipes based on fitting formulae to numerical simula- 
tions or more phenomenological approaches. It yields a rea- 
sonable agreement with numerical simulations and provides 
a competitive approach, since it remains difficult and time- 
consuming to describe a range of scales that spans three 
orders of magnitude or more by ray-tracing simulations. 
Even though this particular model may still be improved, 
these results already show that building on systematic and 



physically motivated models of the 3D matter distribution 
is a promising route to predict both 3D and 2D (i.e. pro- 
jected) statistics. 

An advantage of this approach, as compared with nu- 
merical simulations or global fitting formulas, is that it pro- 
vides a clear link between the observed weak-lensing quanti- 
ties, such as the convergence power spectrum or bispectrum 
on a given range of multipoles £, and the various ingredi- 
ents that govern the underlying 3D matter density field. In 
particular, we can distinguish the scales associated with the 
perturbative regime from those that probe nonperturbative 
features described for instance by the 1-halo term. This is 
useful since the perturbative regime can be predicted by 
rigorous and systematic approaches (e.g., in this paper we 
use both standard perturbation theory and a peculiar re- 
summation scheme, up to one- loop order), which offers a 
well-controlled constraint on the cosmological parameters. 
The nonperturbative regime relies on a more phenomeno- 
logical approach, i.e. a halo model, so that the link to cos- 
mological parameters is somewhat weaker. Indeed, the ac- 
curacy is lower (because there is no rigorous and system- 
atic framework) and the model involves some intermediate 
parameters, such as halo profiles, that show a weak depen- 
dence on cosmology. Therefore, it can be useful to separate 
both regimes when deriving constraints on cosmology from 
observations. Moreover, the explicit link provided by such 
models allows one to use gravitational lensing effects to 
probe the halo properties themselves. 

We find that going to one-loop order brings a more 
significant improvement, with respect to lowest-order per- 
turbation theory, for the lensing bispectrum than for the 
power spectrum. Moreover, while for the bispectrum the 
combination of the 1-halo and 2-halo terms with the one- 
loop perturbative contribution (associated with the 3-halo 
term) provides a reasonable match to simulations, for the 
power spectrum there remains an intermediate regime that 
is not well described by the simplest combination. This sug- 
gests that higher orders of the perturbative expansion play 
a greater role for the power spectrum and should be taken 
into account, or that the matching between the 2-halo and 
1-halo contributions need be improved. In our case, we used 
a simple geometric interpolation for the 3D power spectrum 
and we find that it also provides good results for the weak 
lensing power spectrum. These features agree with previous 
results for the 3D density power spectrum and bisp ectrum 
([Sefusatti et al.ll2010t [Valageas fc Nishimichill2011bD . 

As in 3D, the 2-halo term does not significantly con- 
tribute to the lensing bispectrum for equilateral configura- 
tions. It only plays a role for very "squeezed" configura- 
tions on weakly nonlinear scales, with a ratio of 100 or 
more between the two long sides and the smallest side of 
the Fourier-space triangle. This is also the regime where 
we find a discrepancy between our predictions and the 
simulations. Since both the 3-halo and the 1-halo terms 
provide good agreement with the simulations for all other 
regimes, where they dominate, this suggests that this dis- 
crepancy arises from a lower accuracy of the 2-halo term. 
This is not really surprising since this contribution is a pri- 
ori more intricate, because it combines internal halo proper- 
ties with larger-scale correlations. Nevertheless, since this 
regime corresponds to extreme configurations it is not a 
very serious practical problem. 

Overall, we find that our model provides a good agree- 
ment with the numerical simulations for a variety of cos- 
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mologies and gives a robust framework. It could still be 
improved in various manners. First, the accuracy of the per- 
turbative contribution may be increased by including higher 
orders beyond one-loop or by using alternative resumma- 
tion schemes. Second, the under lying halo model could be 



refined to include substructures (lShethll2003t iGiocoli et al.l 



20101) ■ deviations f rom spherical profiles ( Jing fc Sutoll2002l: 



Smith et al.ll2006f ). or the effect of baryons ( Guillet et al.l 
201011 . Next, the model could be generalized to non- 
Gaussian initial conditions, which should yield distinc - 
tive signatures in the bispectrum (jSefusatti et al.l |2010D . 
Another interesting generalization would be to include the 
effects of massive neutrinos. Indeed, it is still difficult to 
model the damping associated with neutrino free-streaming 
using numerical simulations and simple fitting formulae, 
and analytic or semi- analytic prescriptions cq t ild be use- 
ful (iLesgourgues et al.ll2009t ISaito et al.ll2009t iBird et al.l 
I20m 

In this first study we focused on Fourier-space statis- 
tics, the convergence power spectrum and bispectrum, from 
which other statistics such as real-space correlation func- 
tions can be derived. We shall present our results for such 
real-s pace statistics in a companion paper (jValageas et al.l 
[2MI . 
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Appendix A: Six alternative cosmologies 



We give in Table ETD the six alternative sets of cosmological 
parameters that we consider in Sect. [7] to test the robust- 
ness of our model. The first line is our fiducial cosmology, 
used in other sections, while lines two to seven give the 
cosmological parameters associated with the six scenarios 
where rt^, Q,ch?, and wq are modified by ±10%. The pa- 
rameters marked as "-" remain equal to the fiducial values. 
Several parameters in Table lA.ll may simultaneously vary 
because we use the following set of independent parameters: 
As, Us, Oc/i^, f^b/i^, ride, Wo, and we assume the Universe is 
flat. For instance, if we change flch'^ keeping other param- 
eters fixed, h has to be varied from the following equation; 
h = ^n^K^/n^ = y(fV?TTW)7(l~~^Uj- Then, 
and rib have to be varied as well, because of fib — ^hh^ /h^ 
and the assumption of flat Universe. In the case of varying 
Us and Wo, is changed because other parameters (includ- 
ing As) are fixed. 
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h 




C8 


Wo 


fiducial 


0.238306 


0.761694 


0.0415995 


0.732 


0.958 


0.759208 


-1 


ns+ 










1.0538 


0.920651 




Us — 










0.8622 


0.627354 










0.038429 


0.7616 




0.8149102 










0.045347 


0.7011 




0.7014909 




wo+ 












0.7813410 


-1.1 


Wo— 












0.7329135 


-0.9 



Table A.l. The cosmological parameter values for fiducial and six alternative cosmologies. The parameters marked as 
"-" remain equal to the fiducial values. 
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